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Abstract. We have compiled optical and radio astrometric data of the microquasar LS 5039 and derived its 
proper motion. This, together with the distance and radial velocity of the system, allows us to state that this 
source is escaping from its own regional standard of rest, with a total systemic velocity of about 150 km s~^ 
and a component perpendicular to the galactic plane larger than 100 km s~^. This is probably the result of an 
acceleration obtained during the supernova event that created the compact object in this binary system. We have 
computed the trajectory of LS 5039 in the past, and searched for OB associations and supernova remnants in 
its path. In particular, we have studied the possible association between LS 5039 and the supernova remnant 
G016.8— 01.1, which, despite our efforts, remains dubious. We have also discovered and studied an Hi cavity in 
the ISM, which could have been created by the stellar wind of LS 5039 or by the progenitor of the compact object 
in the system. Finally, in the symmetric supernova explosion scenario, we estimate that at least 17 Mq were lost 
in order to produce the high eccentricity observed. Such a mass loss could also explain the observed runaway 
velocity of the microquasar. 
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1. Introduction 

Microquasars are stellar-mass black holes or neutron stars 
that mimic, on smaller scales, many of the phenomena 
seen in AGN and quasars. These objects have been found 
in X-ray binary systems, where a compact object ac- 
cretes matter from a companion star. Radio emitting 
X-ray binaries with relativistic radio jets, like SS 433, 
GRS 1915-H105, GRO J1655-40 or Cygnus X-3, are good 
examples of microquasars (see Mirabel & Rodriguez 1999 
for a detailed review). With the recent addition of LS 5039 
(Paredes et al. ^OOOD , Cygnus X-1 (Stirling et al. |200l[ ) and 
XTE J1550-564 (Hannikaincn et al. |200l[ ) to the micro- 
quasar group, the current number of this kind of sources 
is 14, among ~ 280 known X-ray binaries (Liu et al. 2000; 
Liu et al. |200l| ) of which ~ 50 display radio emission. 
Recent studies of microquasars can be found in Castro- 
Tirado et al. ( |200l| ). 

Attention to the star LS 5039 was first called by 
Motch et al. ( |1997D , who proposed it as a High Mass X- 
ray Binary (HMXB) candidate associated with the X-ray 
source RX J1826.2— 1450. The object is located at an es- 
timated distance of ^ 3.1 kpc and close to the galactic 



plane {I = 16.88°, b = -1.29°). Soon after, non-thermal 
and moderately variable radio emission was reported by 



Marti et al. (1998) using the Very Large Array (VLA). 



The evidence of its microquasar nature was provided by 



Paredes et al. (2000) when radio jets were discovered with 



the Very Long Baseline Array (VLBA) at milliarcsecond 
(mas) angular scales. These authors also pointed out the 
possible connection of LS 5039 with 3EG J1824-1514, 
i.e., one of the unidentified EGRET sources of high energy 
7-rays . X-ra y observations of RX J1826.2— 1450 by Ribo 
et al. (1999) did not reveal pulsations and were consistent 
with a significantly hard X-ray spectrum up to 30 keV, 
with a strong Gaussian iron line at 6.6 keV. 

The mass donor in LS 5039 was originally classified as 



an 07V((f)) star by Motch et al. ( |1997D . This classification 
has been recently improved thanks to optical and near 
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infrared spectroscopic observations by Clark et al. (2001), 
which indicate an 06.5V((f)) star. 

Our present knowledge of the system orbit is based 
only on the radial veloc ity measurements in the optical 
by McSwain et al. (2001). Their most remarkable findings 
consist of a short orbital period of P — 4.117 days and a 
significant eccentricity of e = 0.41. They were also able to 
determine the system radial velocity and a mass function 
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Table 1. Compilation of optical and radio positions, with associated errors, of LS 5039. 



Wavelength Epoch a (ICRS) cTacosi S (ICRS) ag Catalog 

Domain (h, m, s) (mas) (°, ', ") (mas) name 



Optical 



1905.45 

1907.43 

1951.577 

1979.484 

1986.653 

1991.75 

2000.289 



18 26 
18 26 
18 26 
18 26 
18 26 
18 26 
18 26 



15.0194 

15.0177 

15.034 

15.0557 

15.054 

15.0427 

15.0563 



255 
255 
250 

64 
300 
149 

13 



-14 50 53.300 
-14 50 53.075 
-14 50 53.59 
-14 50 54.075 
-14 50 54.29 
-14 50 54.229 
-14 50 54.277 



247 
247 
250 

47 
300 
120 

13 



AC 2000.2 
AC 2000.2 
USNO-A2.0 
TAG 2.0 
CSC 1.2 
Tycho-2 
UCACl 



Reference 



Urban, private communication 
Urban, priva te com munication 
Monet et al. h999 



Zacharias & Zacharias 
Morrison et al. 



H0g et al. 200C 



2001 



1999 



Zacharias et al. 200C 



Radio 



1998.24 
2000.42 



18 26 15.056 
18 26 15.0566 



10 

4 



-14 50 54.24 
-14 50 54.261 



10 
6 



VLA obs. 
VLBA obs. 



Marti et al [l998 
Ribo et al. EOOT 



of /(m) = 0.00103 Mq. The fact that the available optical 
photometry does not show an ellipsoidal modulation could 
indicate that the orbital inclination is very low. Therefore, 
the possibility of having a black hole in the system cannot 
be ruled out in spite of the small mass function observed. 

The formation of the compact object in an X-ray bi- 
nary necessarily requires a supernova explosion that does 
not disrupt the system. This explosive event is expected 
to considerably change the kinematic properties of the bi- 
nary star. The change may be rather extreme (i.e., kick 
velocities approaching ^ 10'^ km s^^) for highly asym- 
metric supernovae, as it has been proposed for the mi- 
croquasar Circinus X-1 (Tauris et al. 199S| ). This mech- 
anism may also be responsible for ejecting X-ray bina- 
ries into the halo of the Galaxy. The fast moving X-ray 
nova XTE Jl 1184-480 is a possible example (Mirabel et al. 



2001), although it could have been ejected from a globular 
cluster in the past. On the other hand, the microquasars 
GRO J1655— 40 and Cygnus X-1, also disp lay ru naway 
velocities (Shahbaz et al. |l999| ; Kaper et al. |l999|) . With 



independence of their origin in the galactic plane or in the 
halo, the existence of runaway microquasars is a new and 
important issue that deserves an in-depth study. 

In this paper we focus on the kinematic properties 
of LS 5039 and its surroundings. In Sect. |^ we estimate 
the proper motions of the system, using optical and radio 
data. In Sect. ^ we discuss the distance and space velocity 
of LS 5039, whereas in Sect. ^ we analyze the trajectory 
of the binary system in the past. In Sect. || we study the 
possible association of this microquasar with a supernova 
remnant, and in Sect. ^ we analyze their Hi surround- 
ings. Finally, we make a global discussion in Sect. and 
summarize our main conclusions in Sect. p. 



2. Positions and proper motions 

2.1. Optical positions 

LS 5039 {V ~ 11.2, B ~ 12.2) is not included in most ac- 
curate astrometric catalogs, like Hipparcos and Tycho-1 
(which are complete up to F ~ 9.0 and V ~ 10.0, respec- 
tively). Although it appears in Tycho-2, its position has 



a large uncertainty when compared to the average error 
within its magnitude range (H0g et al. 2000). As a result, 
the catalogued proper motions lack of the required pre- 
cision. Moreover, other astrometric catalogs which were 
not included in Tycho-2 estimate of the proper motions, 
have been released since then. In view of all these facts, we 
decided to carry out a thorough study of all available cat- 
alogs. First, we did a search using the query forms within 
the VizieR database (Ochscnbein et al. 2000|) at the Centre 
de Donnees astronomiques de Strasbourg (CDS), and af- 
terwards we looked for new information on several public 
web sites. The results of this search are listed in Table 0, 
with all coordinates in ICRS, and commented below. 

The oldest astrometric information for LS 5039 is 
dated around 1905 and comes from the Astrographic 
Catalog 2000 (AC 2000, Urban et al. |l998D. However, a 



re-reduction of this catalog was performed recently, us- 
ing an improved reference catalog that allowed a better 
handling of the systematic errors (S. E. Urban, private 
communication). This new catalog will soon be released 
as AC 2000.2. Since changes were most significant in the 
faint stars of the southern AC 2000 zones, the position of 
LS 5039 suffered a noticeable shift. Actually, two positions 
corresponding to different epochs are given in AC 2000.2. 

The following available astrometry, dated around 1950, 
comes from the scanning and reduction of the Palomar 
Observatory Sky Survey plates, carried out by Monet et al. 
( |1999D , which is part of the USNO-A2.0 catalog. 

We inspected the Second Cape Photographic Catalog 
(CPC2, Zacharias et al. 1999| ), a Southern Hemisphere as- 
trometric catalog, containing observations between 1962 
and 1973. Unfortunately, LS 5039 is not listed in it, be- 
cause the limiting magnitude of this catalog was V ~ 10.5. 

The Twin Astrographic Catalog version 2 (TAG 2.0) 
provides positional information around 1980. Apart from 
the internal error of CTqcosA' = 62 mas and as — 45 mas 
associated to the LS 5039 position, an external error of 
15 mas has been taken into account following an estimate 
carried out by Zacharias & Zacharias ( 1999 ). 

The last photographic catalog used is the Guide Star 
Catalog version 1.2 (CSC 1.2, Morrison et al. |200l|) . This 
is a re-reduction of the GSC (Lasker at al. 1990 ), after re- 
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moving plate-based systematic distorsions that are a func- 
tion of magnitude and radial distance from the plate cen- 
ter. In addition, by using the Starlink library SLALIB, we 
have transformed the catalogued position from FK5 to 
ICRS. 

As mentioned above, L S 503 9 appears as well in the 
Tycho-2 catalog (Il0g et al. 2000 ). The average internal er- 
ror in position for a V ~ 11.2 star is around 50 mas, while 
the errors in a cos 5 and 5 for LS 5039 are 93 and 120 mas, 
respectively. This fact tells us that the position estimate 
is not as good as one would expect. On the other hand, 
the scatter-based errors from the 4 positions present in the 
Tycho-2 data are 149 and 60 mas, respectively. Hence, we 
have considered CTqcosiS = 149 mas and as = 120 mas. 

A second epoch of the GSC was performed during the 
nineties. However, the positional errors in the current ver- 
sion of this catalog, GSC 2.2.01, are intended as indicators 
for operational use only, and cannot be used for scientific 
studies. Hence, we have considered that this catalog does 
not provide accurate astrometric information to be used 
for estimating proper motions. 

Finally, we have used the first release of the US 
Naval Observat ory C CD Astrograph Catalog (UCACl, 
Zacharias et al. 200C). As can be seen in Table |l|, this 
is by far the best optical astrometric catalog available. 

2.2. Radio positions 

Independent astrometric information can be obtained 
from radio observations of LS 5039. The first available as- 
trometric information comes from the NVSS by Condon 
et al. ( 199§| ). However, the position error in this catalog, 
around 1", is too large to be of any use for our purposes. 
Martf et al. ( |1998 ) carried out VLA^A configuration ob- 
servations which provided accurate (10 mas) astrometry 
for LS 5039. Finally, phase-referencing VLBA^ observa- 
tions conducted by Ribo et al. ( 2002D provide the last ra- 
dio position, with the best accuracy among all available 
data. Details of the astrometry for the last two cases are 
listed in Table ^ The relatively large location error of few 
mas in the VLBA data is due to confusion produced by 
the high level of Galactic electron scattering in this region 
(for both the reference source and LS 5039). 

We must state that no correction for parallax has been 
applied to the obtained positions, since at an estimated 
distance of ~ 3 kpc, this would translate into less than half 
a mas, which is always much smaller than the available 
uncertainties. 



2.3. Proper motions 

It is clear from the astrometric data that both, the optical 
and radio sources, are almost in the same position of the 



Table 2. Proper motions estimates for LS 5039. 



^ The VLA and the VLBA are operated by the National 
Radio Astronomy Observatory (NRAO). The NRAO is a facil- 
ity of the National Science Foundation operated under coop- 
erative agreement by Associated Universities, Inc. 



Data set 


f-^a cos 5 






(mas yr~^) 


(mas yr~^) 


Optical 


4.7 ± 1.3 


-11.0 ±0.8 


Radio 


4.0 ±4.9 


-9.6 ±5.3 


Optical-f Radio 


4.7 ± 1.1 


-10.6 ± 1.0 



sky. However, in order to show that this is not a chance 
coincidence and that both emissions originate in the same 
object, we have computed independent proper motions for 
the optical and radio data. 

Although the position uncertainties from the 
AC 2000.2 and USNO-A2.0 catalogs are relatively 
large, the epoch span achieved justifies their inclusion 
in the estimate of the proper motions. Hence, from 
the optical data, and taking into account the astro- 
metric uncertainties, we obtain the following results: 
Ma cos(5 = 4.7 ±1.3 mas yr~^, = — 11.0 ± 0.8 mas yr~^, 
where the errors come directly from the least squares fit. 
The two accurate radio positions give proper motions of: 
A*Qcos(5 — 4.0 ± 4.9 mas yr^^, — —9.6 ± 5.3 mas yr~^. 
Although this last result has a large error because 
only two points are available, we can say that both, 
the optical and radio sources, have very similar proper 
motions. Therefore, based only on astrometric data, 
we are able to confirm that both, the optical and the 
radio emission, originate in the same object. Hence, 
we will use all the data listed in Table |l] to compute 
accurate proper motions for LS 5039, which happen to be: 
A*Qcos(5 — 4.7± 1.1 mas yr~^, — — 10.6± 1.0 mas yr~^. 
All these results are summarized in Table |^. We can 
now transform the proper motions into galactic co- 
ordinates and obtain fii = —7.2 ± 1.0 mas yr^^ and 
/if, = —9.1 ± 1.0 mas yr^^. It is clear from these results 
that there is a noticeable motion perpendicular to the 
galactic plane and moving away from it. 

According to our least squares fits, and defining t = 
yr — 2000.0, the predicted ICRS values for a and S near 
t = are: 



18'' 26" 15.0565" 



-14° 50' 54.260" 



A.7t/ cos(5 



mas 



10.6 1 



mas 



(1) 



(2) 



being 1^9 + 1.2 1'^/ cos^ S and \/9 + t'^ the errors in mas 
in a and (5, respectively. Comparing the fits with the po- 
sitions in Table ^ and the proper motions in Table ^ we 
can say that, approximately, the offsets are determined by 
the VLBA data point because of the small error in posi- 
tion, while the proper motions are obtained by the optical 
points due to the huge time span. 

In Fig. |l| we have plotted the offsets in a cos S and S 
from the fitted values for year 2000.0 versus time, for the 
data listed in Table |l|, together with the respective fits to 
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Fig. 1. Offsets in Right Ascension (left) and Declination (right), relatives to the fitted ones for year 2000.0, versus 
time for all the positions listed in Table |f|. Open circles represent the optical positions, and filled circles the radio ones. 
The solid lines represent the least squares fits to the whole data sets. 



the data. Notice that in both plots the UCACl and VLB A 
data are almost superimposed. 



3. Ejection from the galactic plane 

Since the distance to LS 5039 and its uncertainty are fun- 
damental to compute the space velocity of the system, we 
think it is worth to perform first of all a short study on 
this issue. 



3.1. Distance and its uncertainty 



Motch et al. (1997) proposed a distance of 3.1 kpc to 
LS 5039 based on the star color excess, and stated that this 
estimate could have a large uncertainty. This value was 
computed assuming an 07V((f)) spectral type. However, 
recent observations by Clark et al. (2001) show that the 
optical companion has a spectral type of 06.5V((f)). On 
the other hand, old calibrations for the intrinsic color in- 
dex a nd ab solute magnitude (Johnson 1966 , Deutschman 

(|1997|) . Hence 



we 



et al. |l976| ) were used by Motch et al 
have performed a new estimate of the distance taking into 
account the new spectral type and more recent calibra- 
tions. 

The optical photometry available up to now, contain- 
ing at least B and V magnitudes comes from Drilling 
( |l99l| ), LahuUa & Hilton ( |l99^ ) and Clark et al. ( pOOl| ), 
and is listed in Table ^. Using an intrinsic color index of 
{B - V)a ^ -0.30 ± 0.02 for an 06. SV star (Schaerer 
et al. |l99e| , Lejeune & Schaerer |200lD , we can com- 
pute the color excess E{B — V) for all the observations. 
Finally, using the relationship Ay = (3.30 + 0.28 (^ - 
V)o + 0.04 Eb-v)E{B - V) (Schmidt-Kaler |l982t) and 
My = -4.99 ± 0.3 for an 06.5V star (Vacca et al . (l996|) 



Table 3. Distance estimates to LS 5039 derived from op- 
tical photometry. 



V 


B-V 


d (kpc) 


1L23±0.01'' 


0.95 ±0.01" 


2.68 ± 0.06 


11.20 ±0.02'' 


0.95 ±0.02'' 


2.64 ±0.06 


11.33 ±0.02'= 


0.85 ±0.02'= 


3.26 ± 0.07 


11.32 ±0.01'' 


0.85 ±0.02'' 


3.25 ± 0.07 



Drilling (1991) 



" LahuUa & Hilton ( |1992| ) 

'= 1996 October observations by Clark et al. ( |200l|) 
'' 2000 September observations by Clark et al. ( |200lD 



we obtain the distance estimates listed in Table y. The 
weighted mean of these values is 2.9 ± 0.3 kpc, and will 



be the distance to LS 5039 used hereafter. Finally, assum- 
ing that the Sun is at 8.5 kpc from the galactic center, 
we obtain a galactocentric distance of 5.8 ± 0.3 kpc for 
LS 5039. 



3.2. Space velocity 

If we assume a distance of 2.9 ± 0.3 kpc to LS 5039, 
as seen in the previous section, and a systemic VJ- = 
4.6 ± 0.5 km (McSwain et al. |2001| ), the proper mo- 
tions estimates translate into {U — 40±5, V = — 82± 16, 
W = -118 ± 19) km s^i in the Local Standard of 
Rest (LSR) defined by (Uq = 9, Vq = 12, Wq = 7) 
km s^^. This gives a total systemic velocity of Vgys — 

of 215 km s~^ for the galactic rotation at 5.8 kpc from 
the galactic center (Fich et al. 1989) we can transform 
the galactic plane space velocities {U and V) from the 
LSR into the LS 5039 regional standard of rest (RSR). 
Applying this transformation we find U = 51 ± 6 and 



M. Ribo et al.: LS 5039: a runaway microquasar ejected from the galactic plane 



5 



Table 4. Space velocity estimates, in km s ^, for LS 5039 
related to the LSR and to its RSR. 



Frame 


U 


V 


W 




LSR 


40 ±5 


-82± 16 


-118 ± 19 


149 ± 18 


RSR 


51 ±6 


-71 ±16 


-118 ± 19 


147 ± 17 



V = -71 ± 16 km s-\ and Vsys = 147 ± 17 km s^^. All 
these results are summarized in Table ^. 

Taking into account that the cosmic dispersion for 
early-type stars is (ct(7, cry, o-vf) — (7,8,4) km s^-'^ (Torra 
et al. pOOO| ) we can conclude that this source is escaping 
from its own RSR and has a large velocity component per- 
pendicular to the galactic plane. Since LS 5039 is a HMXB 
containing an early-type star, it seems very unlikely to be a 
high speed halo object crossing the galactic plane. Hence, 
the most plausible explanation for the observed velocity 
is that the binary system obtained an acceleration during 
the supernova event that created the compact object in 
this microquasar. 



4. The past trajectory of LS 5039 

An interesting thing to do once a position and space ve- 
locity estimates are available, is to compute the trajectory 
of LS 5039 in the past, and then to look for OB associ- 
ations and SN Remnants (SNRs) in its path, in order to 
establish possible relationships. For this purpose, we com- 
puted the galactic orbital motion of this system under the 
gravitational field of the Gala xy. The galactic mass model 
by Dauphole & Colin (1995) was adopted for the inte- 
gration, which was was performed using a Runge-Kutta 
fourth-order integrator using time steps of 1000 years. 

It is helpful to search for OB associations in the past 
trajec tory o f runaway X-ray binaries, as done by Ankay 
et al. ( |2001| ) with HD 153919/4U 1700-37, because if ev- 
idences are found that the runaway system originated in 
an OB association, it is possible to estimate the age of the 
binary system after the SN explosion. According to the 
computed trajectory in the past for LS 5039, there are 
two OB associations in or close to its path in the plane 
of th e sky: Set 0B3 and Ser 0B2 (Melnik & Efremov 
1997). Unfortunately, the distance to Set 0B3 is ^ 1.5 kpc 
(Melnik & Efremov |1997|), a nd the distance to Ser 0B2 is 
1.9 ± 0.3 kpc (Forbes 2000). Hence, the two OB associa- 
tions found in the path of LS 5039 are too close to us to 
be related to it. 

Since we have not been able to fix a limit on the 
integration time in the past thanks to a possible rela- 
tionship with an OB association, we will use the ver- 
tical distribution of early-type stars in the Galaxy for 
this purpose. 0-type stars are typically located at dis- 
tances within 45 ± 20 pc from the galactic midplane (Reed 
2000| ). According to the galactic latitude of LS 5039 and 
using a distance of 2.9 ± 0.3 kpc, its actual height is 
Z = —65 ± 7 pc. Since it is escaping from the galactic 



plane, we can compute the trajectory backward in time up 
to when it had a height of Z = 65 pc, which is ^ 1.1 Myr 
ago. Hence, the SN explosion probably took place in the 
last 1.1 Myr. 

Thus, for a possible association of LS 5039 with a SNR, 
we will focus our attention on the galactic trajectory of 
the binary system for the last ~ 1 Myr. In addition, the 
likelihood of detecting a SNR decreases with time, so it is 
not expected that the radio emission of the SNR can be 
observed much beyond the time interval considered above 



(e.g., Shklovskii|l96§). 

A s earch in a catalog of galactic supernova remnants 
(Green 200C ) reveals the presence of several SNRs near the 
path of LS 5039 on the plane of the sky, though the center 
of none of them is crossed by its trajectory. Only three of 
them, namely SNR G016.7-f 00.1, SNR G16. 8-01.1, and 
SNR 0018.84-00.3, are within 1° of the LS 5039 path. 
The distance to SNR G016.7-f 00.1 was estimated to be at 
~ 14 kpc by Reynoso & Manguni ( ^000 ), maki ng imp ossi- 
ble an association with LS 5039. Dubner et al. (1999) con- 
cluded that SNR G018.8-h00.3 is located at 1.9 ±0.5 kpc, 
and has an age of ~ 16 000 yr. Although the distance is not 
very different from the one to LS 5039, the age of this SNR 
is by far too short to be associated with the microquasar. 
Hence, the only remaining candidate is SNR G16.8— 01.1, 
which deserves a detailed study. 

5. SNR G016.8-01.1 

A wide field radio map of the surrounding of LS 5039 is 
shown in Fig. ^ with the microquasar (contours) and the 
nearby SNR G016.8— 01.1 (grey scale), together with the 
position of the Hll-region RCW 164 (which is the stronger 
source in the field). The arrow in this figure marks the 
proper motion sense of LS 5039 as it would be seen from 
the LSR (i.e., correction for the peculiar velocity of the 
Sun has been applied), while the dashed line represents 
the trajectory in the past up to 10^ yr, with the corre- 
sponding error boxes in position at that epoch and 5 10** yr 
ago. From this figure we can see that LS 5039 crosses the 
projection in the plane of the sky of SNR G016.8— 01.1 
between ~ 4 10"* and ^ 1.3 10^ years ago. 

If one of the two stars forming a close binary system 
experiences a SN explosion, it may form a SNR and an 
X-ray binary. Since conservation of the linear momentum 
after the SN explosion is required, if the SNR and the X- 
ray binary are related, we expect them to be aligned with 
the proper motion of the latter, and any deviation from 
this behavior should be explained by the projection in the 
plane of the sky of the initial space velocity of the binary 
system prior to the SN event. This does not seem to be the 
case here, because SNR G016.8— 01.1 and LS 5039 are not 
well aligned with the LS 5039 proper motion. Moreover, 
if we assume that the center of the shell-like remnant is 
located at the minimum of radio emission in Fig. ^ ap- 
proximately located at a = 18*^ 25"^ 03", (5 = -14° 37' 30" 
(or / = 16.94°, b = —0.93°), a peculiar velocity in the 
plane of the sky for the original system of ~ 70 km s^^ 
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Fig. 2. Wide field radio map of LS 5039, its nearby 
shell-like SNR G016. 8-01.1 and the H ii-region RCW 164 
(which is the stronger source in the field). The grey scale 
emission is taken from the Parkes-MIT-NRAO tropical 



survey at the 6 cm wavelength (Tasker et al. 1994). The 
overlaid contours correspond to the NVSS map of the re- 
gion at the 20 cm wavelength (Condon et al. |1998 ). The ar- 
row marks the proper motion sense (see text). The dashed 
line is the computed trajectory for the last 10^ yr, with the 
corresponding error boxes in position at that epoch and 
5 10^ yr ago. The dotted line represents positions with a 
galactic latitude of —0.8° for reference purposes. 



is obtained. This velocity is much higher than the typi- 
cal ^ 10 km s^^ found for early-type stars. In order to 
reduce this velocity we can use in our favor the error in 
the proper motion angle, and find that it turns out to be 
~ 50 km s^^, which is still too high. However, the cen- 
ter of the shell-like remnant could be in another position, 
closer to the LS 5039 trajectory, and the derived peculiar 
velocity for the system prior to the SN explosion could fit 
in the typical values of early-type stars. Hence, from the 
kinematical point of view, we are not able to firmly discard 
a possible association between both objects. Therefore, we 
have performed an in-depth study of the remnant. 

SN R GO 16. 8-01.1 was discovered by Reich et al. 
( |1986D using a combination of survey data analysis and 
new multifrequency observations. The extended radio 
source has a complex structure due to the superposition 
of the Hil-region RCW 164 (Rodgers et al. |1960|), which 



is coincident with the peak of the continuum emission (in 
black in Fig. ^) . The diffuse radiation obs erved around the 
peak is strongly polarized (Rodgers et al. 1960 ), indicating 
a synchrotron nature. Polarization decreases towards the 
position of the foreground Hii-region as a result of the 
thermal contribution. Thermal and non-thermal compo- 
nents of the radiation cannot be clearly separated, but the 
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Fig. 3. H166q; line observations of the Hii-region 
RCW 164. The solid line represents a Gaussian fit to the 
data, with its maximum being ai v — 16.5 ± 0.8 km s~^. 

total fiux from the SNR at 5 GHz seems to be '--^ 1 Jy (J. A. 
Combi, private communication). The angular diameter of 
the remnant is ^ 30'; its distance remains unknown. 

5.1. A lower limit of the distance 

In order to determine a lower limit of the distance to 
SNR G16.8— 1.10, we carried out H166a recombination 
line (1424.734 MHz) observations of the foreground Hil- 
region in November 17th 2000. We used a 30-m radiote- 
lescope at the Instituto Argentino de Radioastronomi'a 
(lAR), Villa Ehsa. The receiver is a hehum-cooled HEMT 
amplifier with a 1008-channel autocorrelator at the back 
end. The HPBW at a wavelength of 21 cm is 30' and the 
temperature of the system on the cold sky during the ob- 
servations was about 35 K. The H166q; line was detected 
at a velocity of 16.5 ± 0.8 km s^^ (see Fig. ^ after 1.5 
hours of integration time, with a signal-to-noise ratio of 
^ 4. At a location of I « 16.8°, b w —1.1°, standard galac- 
tic rotation models (Fich et al. 



198C) 



indicate that the 
observed velocity corresponds to a distance d ~ 1.8 kpc. 
Consequently, SNR G16.8— 1.10 should be farther than 
~ 2 kpc, as suggested by its relatively small angular size. 



5.2. Particle density estimates 

If we adopt for SNR G016.8— 01.1 the same distance to 
LS 5039, we can make some estimates of the parameters 
that would characterize the SNR in case of being associ- 
ated with the X-ray binary. Using a size of 30' and a dis- 
tance of 2.9 kpc, the inferred radius would be i? = 12.7 pc 
and it would be, thus, still in t he ad iabatic expansion 
phase. Using the standard Sedov (1959) solutions, we can 
express the particle density of the ambient medium as 



4.4410"«t^£;5ii?j; 



(3) 
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where t is the age of the remnant in years, E^i is the orig- 
inal energy release of the SN explosion in units of 10^^ erg 
and Ri is the remnant radius in units of 10 pc. Since the bi- 
nary system is bound after the SN explosion, its total mass 
has to be larger than the mass of the SNR (van den Heuvel 



197S), which, therefore, should be moving at least with the 



same space velocity in the opposite direction. Hence, the 
maximum age of the remnant if it were related to LS 5039, 
would be around 5 10^ yr. If we use this age and assume 
that E^i = 0.4 (Spitzer 1998), we find an ambient density 
of n ~ 13.5 cm""^. Considering the uncertainty in the dis- 
tance, we can say that 8 < n < 23 cm^'^. Although such 
number densities would not be very unusual towards this 
direction of the inner part of the Galaxy, we have made 
H I observations in order to estimate the typical column 
densities in the region of interest. 

The H I observations were performed during November 
29th and December 1st 2000, with the already mentioned 
JAR telescope. The H i line was observed in a hybrid total- 
power mode with a sampling on a 0.25°-lattice around the 
position of SNR G016.8-01.1, covering a total of 4° x 4°. 
The velocity resolution obtained was 1.05 km s^^, with a 
total coverage of ±450 km s~^. The rms level in brightness 
temperature was ~ 0.1 K. Those channel maps in the 
velocity range from 10 to 50 km s~^ (corresponding to 
distances between 1.5 and 4.5 kpc) were inspected looking 
for clouds or a local minimum that could be associated to 
the SNR. 

The maps show a strong gradient of brightness tem- 
perature towards the galactic plane, but no clear evidence 
for a structure that could be associated with the SNR. 
In Fig. ^, we show the integrated column density map 
for the velocity interval between 26-34 km s~^, where the 
contour labels are in units of lO'^^ cm~^. With the column 
densities observed in the region where SNR G016.8— 01.1 
is located (~ 4 lO^*' cm^^), the ambient density should 
be ~ 5 cm~'^. This is nearly a factor of 3 lower than the 
estimated density from the expected size of the SNR if it 
is at 2.9 kpc. Taking into account the uncertainty in the 
distance, the ambient density determined from the rem- 
nant's size (8 < n < 23 cm~'^) would be slightly larger 
than that estimated from the Hi observations. It could 
be that there is more material under the form of other 
molecular species, or that the original SN energy release 
could have been different from what we have assumed, or, 
finally, that the distance to the SNR is not 2.9 kpc and it 
is not related to LS 5039. The case, consequently, remains 
unsolved. 



6. The H I surroundings of LS 5039 

Notwithstanding the absence of structures that can be 
clearly associated with the SNR in the Hi distribution, 
the channel maps in the interval 10-46 km s~^ (see Fig. ||) 
reveal the existence of a large, semi-open cavity which is 
very similar to bubbles blown out by early-type stars in re- 
gions with steep density gradients (e.g., Benaglia & Cappa 
199£). This phenomenon has also been observed around 
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Fig. 4. Neutral hydrogen column density distribution 
towards SNR G016.8— 01.1, integrated over the veloc- 
ity interval 26-34 km s^^. Contour levels are in units 
of 10^^ cm~^. The circle marks the position of SNR 
G016. 8-01.1, while the star represents LS 5039. 



another HMXB (HD 153919) by BenagHa & Cappa ( |l999D . 
The ambient material is thought to be swept up by the 
strong wind of the massive early-type star in the binary 
system creating a local minimum around the star. Density 
gradients towards the galactic plane in the original matter 
distribution can result in a preferred escape direction for 
the wind, yielding open structures. Even in some cases the 
star appears displaced from the actual minimum in the H i 
distribution. 

The bubble detected around LS 5039 is centered at a 
systemic velocity of Vgys — 33 km s""'^, with an expansion 
velocity Wcxp — 12 km a~^. The galactic rotation model by 



Fich et al. (1989) indicates that it is located at a distance 
of ~ 3.2 kpc, quite consistent with the distance to the 
microquasar. The total mass removed in order to create 
the cavity is ~ 1.5 10"* Mq, implying a kinetic energy ~ 
2 lO''^ erg s~^, which can be considered as a rough estimate 
of the energy deposited by the wind of the star into the 
ISM. Such a value is in accordance with previous estimates 
for sim ilar O f star wind blown bubbles (e.g., Benaglia & 
Cappa 1999). With a radius of ^ 60 pc, the dynamical age 
of the bubble is 3 10^ yr, which should be considered as 
an upper limit for its actual age since it ignores possible 
contributions from the SN explosion and the wind of the 
progenitor of the compact star. 

We have also searched for other contributors to the 
cavity, among stars earlier than B2 that were projected 
over the cavity, and whose distances were in agreement 
with that of LS 5039. It seems unlikely that the Wolf Rayet 
star WR 115 (WN6-I-0B?) could contribute in some way. 
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Fig. 5. Neutral hydrogen channel maps for the velocity interval 10-46 km s~^, towards LS 5039. Each map covers 
4 km s~^, the central velocity is given in each label. A local minimum and an open low density cavity can be clearly 
seen close to the HMXB. The cavity is open towards the direction opposite to the strongest mass concentrations near 
the galactic plane, as expected from a wind-blown structure. 



because its distance is 2 kpc, and is located over the shell 
of the cavity {I = 16.98°, b -1.03°). The other lumi- 
nous stars in the field, LS 5005, LS 5017, LS 5022, LS 5047, 
LS 5048 and HD 169673 can sum up at most 20% of the 
energy required to create such a minimum in the Hi dis- 
tribution. The 0((f)) star in LS 5039, then, seems to be 
the main agent forming the bubble. 

However, the actual space velocity indicates that the 
system has been close to the cavity only during the last 
2 10^ yr, which is much shorter than the upper limit 



found for the dynamical age of the bubble (~ 3 10^ yr). 
Nevertheless, the cavity could be originated by LS 5039 if 
the pre-SN system was close to the actual position, in such 
a way that the progenitor of the compact object could have 
also contributed to the formation of the bubble. In such 
a case, the X-ray binary would be only about ~ 2 10^ yr 
old, and we would expect to find the radio SNR origi- 
nated in the SN event with a current surface brightness 
of S408 - 3 10-22 W m-2 Hz-i sr-i (Caswell & Lerche 



1979). Filtering techniques of the diffuse galactic emis- 
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sion (e.g., Combi ct al. 1998) allow to detect very low 
brightness SNRs with S408 10^22 ^ ^-2 ^^-1 ^^-i^ 
but nothing has been found in this region except for SNR 
G016.8— 01.1, whose surface brightness at 408 MHz would 
be E408 ~ 6 10~^^ W m^^ Hz^^ sr^^ for a typical spectral 
index a = —0.5. This level of flux is expected at this lati- 
tude for a younger SNR, say with t ~ 10^ yr, but in such 
a case the particle density of the medium around SNR 
G016.8— 01.1 should be higher than the estimates given in 
Sect. 5.2, with values of n '--^ 50 cm~'^, in order to con- 
fine the SNR to its inferred size at d '--^ 2.9 kpc. These 
densities are 1 order of magnitude higher than those de- 
rived from the H I observations. However, a clear rejection 
of SNR G016.8— 01.1 cannot be established on this basis 
alone, due to the flux contamination from the H 11 region, 
which could be responsible for errors as large as 100 % 
in the flux density estimate of the source. Consequently, 
a picture where the Hi bubble has been created by both 
the stellar winds from the binary plus the SN explosion 
remains as an open possibility. 

7. Discussion 

There are two ways to accelerate a binary system by a su- 
pernova explosion: ejection of material from the binary in 
a symmetric SN or an additional velocity kick, produced 
by asymmetries in the SN itself. In Nelemans et al. ( 1999| ), 
the authors conclude that there is no need of additional ve- 
locity kicks in systems like Cygnus X-1 or GRO J1655— 40, 
both containing black hole candidates as the compact ob- 
jects, to explain their observed large runaway velocities. 
On the other hand, from the study of Be/X-r ay bin aries 
containing neutron stars, van den Heuvel et al. (|2000|) con- 
clude that, in order to explain the observed large eccen- 
tricities and low space velocities of these systems, a kick in 
an asymmetric SN explosion scenario is needed. Moreover, 
the microquasar Circinus X-1, containing a neutron star, 
seems to have experienced a highly asymmetric SN ex- 
plosion (Tauris et al. 199£). Unfortunately, the nature of 



the compact object in LS 5039 still remains unknown. 
However, we can try to study it in the symmetric SN sce- 
nario. 

Let Mx and Mq be the masses of the compact ob- 
ject and the optical companion, respectively. Taking into 
account the stellar evolution tracks by Schaerer et al. 
( |1996D , we can assume that Mq ~ 40 Mq . Considering 
that it must rotate slower than the breakup speed, and 
using the mass fun ction /(m) = 0.00103 ± 0.00020 Mq 
by McSwain et al. ( ^OOl] ), an upper limit of Mx < 9 Mq 
is obtained. The minimum mass for the compact object 
would be that of a neutron star of the Chandrasekar mass, 
Mx > 1.4 Mq. Hence, a total mass of 41-49 Mq re- 
mains in the system. 

Since the actual eccentricity of the system is e = 
0.41 ± 0.05, and has a short orbital period of P = 
4.117 ± 0.011 days, both parameters are expected to de- 
crease with time due to tidal forces which act to re- 
circularize the orbit. Hence, the eccentricity just after the 



SN explosion was Cpost-SN > 0.41. Using the equation 
AM = epost-SN (Mx + Mo) and Mx = 1.4 Mq we find 
that at least 17Mq where lost in the SN explosion in the 
symmetric case. If we adopt Mx = 9 Mq the minimum 
mass loss is 20 Mq. Moreover, in order to have the system 
bound after the explosion, i.e., Cpost-SN < 1, the mass loss 
during the SN event must be less than the remaining mass 
in the system, i.e., < 49 Mq. Hence, the values for the re- 
duced mass of the system, defined as ^ — l/(l + epost-SN), 
range from 0.5 to 0.7. Using Eq. (5) of Nelemans et al. 
( 1999D with the actual values of P and e we find that the 
re-circularized period will be Prc-circ — 3 days. Since the 
period of the initial, pre-SN, binary system is given by 
Pi ~ /i^Pro-circ, the orbital period before the SN explo- 
sion was in the range 0.8 to 1.5 days. Using Kepler's third 
law and assuming the latter value to allow the maximum 
separation possible between the progenitor stars, initial 
total masses between 60-100 Mq lead to semimajor axis 
in the range 22-26 Rq. This means that the two stars did 
not evolve separately, but probably through a common 
envelope phase. Hence, a detailed evolution of the system 
should be carried out in order to try to compute the ini- 
tial mass of the progenitor star of the compact object, 
the maximum age of the system before the SN event, and 
the mass loss during the SN explosion itself. This kind of 
study is beyond the scope of this paper. 

Nevertheless, we can try to explain the measured Wgys 
in the context of the symmetric SN explosion scenario us- 
ing the following equation from Nelemans et al. (1999): 



= 213 



A 



Mc^JWq^ 



Mo\ /P 



day 



Mx+M q 

Mq 



(4) 



= 3 days we 



Assuming that Mq = 40 Mq and 
can find two minimum values of Vgys depending on the 
adopted mass for the compact object. For Mx = 1.4 Mq, 
and hence AM > 17 Mq, we find u^ys Z 200 km s"^ 
Similarly, for Mx = 9Mq, and hence AM > 20 Mq, we 
find Vsys ^ 180 km s~^. Our measured recoil velocity is 



Vsys — 150 ± 20 km s^^, which would be perfectly reached 
in the SN explosion scenario, as we have just seen. In fact, 
the measured value is a little bit lower than the computed 
ones, which could indicate that part of the eccentricity of 
the system was reached thanks to a kick in an asymmetric 
SN explosion. However, we must be cautious with the use 
of the value for Vgys obtained after Eq. ^, because some 
effects have been ignored, as pointed out by Nelemans 
et al. (1999), and also the errors in e and Pro-circ could 
reduce the computed values. Hence, we can say that both, 
the high eccentricity measured by McSwain et al. (2001), 
and the high space velocity reported here, are perfectly 
compatible, within errors, with a symmetric SN explosion 
scenario. In other words, the actual eccentricity of the sys- 
tem is naturally explained with the mass loss necessary to 
explain the measured space (recoil) velocity of LS 5039. 

Another common mechanism referred to produce run- 
away stars is based on dynamical ejection from young open 
clusters. Numerical simulations (Kiseleva et al. 1998) show 
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that it can explain some moderately high- velocity stars ob- 
served in the Galactic disk, but only about 1% of the stars 
would be ejected with velocities higher than 30 km s^^. 
Hence, it seems difficult to explain the high observed ve- 
locity of LS 5039 as a result of such a mechanism. 

It is interesting to note that, contrary to other HMXB 
with low orbital periods, LS 5039 is still in the circulariza- 
tion process. Using the actual value of P, a typical value of 
R = 10.3 Rq for the radius of the 06.5V star, and differ- 
ent formalisms (Claret et al. 1995; Claret & Cunha 1997) 



we obtain typical circularization timescales of the order 
of (0.5-5)10^ yr. These values are compatible with the 
system being still in the circularization process, since its 
formation just in the galactic midplane would take 5 10^ yr 
to bring it to the actual position according the integration 
of its trajectory. 

At this point we might ask what is the point of study- 
ing runaway microquasars. In the case of LS 5039, if the 
SN explosion had taken place recently, we could expect 
the system to survive for the following ~ 5 Myr, because 
the 06.5 star is still in the main sequence. Hence, accord- 
ing to the computed trajectory, it would reach a height 
of Z ~ —600 pc, which translates into b ~ —12°. On the 
other hand, LS 5039 could be associated with the high 
energy 7-ray source 3EG J1824— 1514, as suggested by 
Paredes et al. ( 2000| ). Hence, if this association is correct, 
we could be able to detect 7-ray microquasars up to val- 
ues of |6| > 10°, although all confirmed microquasars lie 
within |6| < 5°. In particular, runaway microquasars could 
be connected with some of the unidentified faint, variable, 
and soft 7-ray EGRET sources above the galactic plane , 
as suggested by Romero ( 2001 ) and Mirabel et al. ( 2001 ). 

Hence, it would be interesting to study the proper mo- 
tions and radial velocities of as many microquasars as pos- 
sible. Unfortunately, this kind of information is not easy 
to obtain for a wide sample of systems. Proper motions 
can be obtained either from optical data or either from 
radio data. However, most optical counterparts of these 
sources are faint objects not present in old astrometric 
catalogs. Hence, it is mandatory to acquire new positions. 
In this context, it is better to do it in the radio domain 
because the uncertainties are always smaller than those 
obtained in the optical. On the other hand, it is necessary 
to carry out optical spectroscopy to obtain the radial ve- 
locity curve of these sources and use the correct value for 
the radial velocity of the system. 

An alternative approach could be the search for signa- 
tures of bow shocks in the microquasar vicinity. However, 
recent studies by Huthoff & Kaper ( 2002] ) of runaway OB 
stars, reveal that the success of this method is highly de- 
pendant on the density of the ISM around the runaway 
object. In particular, in only one object, namely Vela X-1, 
a bow shock has been detected. Finally, one could search 
for signatures of explosive events in the ISM. Sensitive 
spectral line observations in the radio are probably the 
best tool for this p urpose, as done in GRO J1655— 40 by 
Combi et al. (|200l|) . 



In any case, the study of runaway microquasars such 
as LS 5039 is likely to contribute significantly to different 
areas of modern high-energy Astrophysics with indepen- 
dence of the observing technique. 

8. Conclusions 

After an in-depth study of the proper motions and sur- 
romidings of LS 5039 our main conclusions are: 

1. Positions at optical and radio wavelengths have been 
used to compute independent optical and radio proper 
motions, which are perfectly compatible. Therefore, 
based only on astrometric data we are able to confirm 
that both, the optical and the radio emission, originate 
in the same object. 

2. From the combined optical and radio positions we 
have computed an accurate proper motion for LS 5039. 
This, together with the new estimate of 2.9 ± 0.3 kpc 
for the distance, allows us to compute a space veloc- 
ity of {U = 51, y = -71, W = -118) km s'^ in its 
Regional Standard of Rest (RSR) . This results implies 
that LS 5039 is a runaway microquasar with Vsys ~ 150 
km s~^, escaping from its own RSR with a large ve- 
locity component perpendicular to the galactic plane. 
This is probably the result of the SN event that created 
the compact object in this binary system. 

3. We have computed the past trajectory of LS 5039. Two 
OB associations have been found close to its path in 
the plane of the sky. However, they are too close to 
us to be related to the microquasar. On the other 
hand, we have also found three SNRs near the path 
of LS 5039. After discarding two of them based on 
distance arguments, we have focused our attention on 
SNR G016. 8-01.1. A study of this source could not 
clearly confirm nor reject the association due to the 
large uncertainties in the estimated radio flux density 
of the SNR. This fact perhaps justifies future, high 
sensitivity searches of low-brightness remnants in this 
region. 

4. We have found a semi-open Hi cavity close to the 
LS 5039 position. Although the 0((f)) star in this mi- 
croquasar seems to be the main agent forming the bub- 
ble, a contribution from the progenitor of the compact 
object cannot be ruled out. 

5. Finally, we are able to explain both, the high space ve- 
locity and the high eccentricity observed, in a symmet- 
ric SN explosion scenario with a mass loss of A M ~ 
17 Mq. 
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